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Abstract. The two-layer quasigeostrophic flow model is an intermidi- 
ate system between the single-layer 2D barotropic flow model and the 
continuously stratified, 3D baroclinic flow model. This model is widely 
used to investigate basic mechanisms in geophysical flows, such as baro- 
clinic effects, the Gulf Stream and subtropical gyres. The wind forcing 
acts only on the top layer. We consider the two-layer quasigeostrophic 
model under stochastic wind forcing. We first transformed this system 
into a coupled system of random partial differential equations and then 
show that the asymptotic probabilistic dynamics of this system depends 
only on the top fluid layer. Namely, in the probability sense and asymp- 
totically, the dynamics of the two-layer quasigeostrophic fluid system is 
determinied by the top fluid layer, or, the bottom fluid layer is slaved 
by the top fluid layer. This conclusion is true provided that the Wiener 
process and the fluid parameters satisfy a certain condition. In partic- 
ular, this latter condition is satisfied when the trace of the covariance 
operator of the Wiener process is controled by a certain upper bound, 
and the Ekman constant r is sufficiently large. Note that the generalized 
time derivative of the Wiener process models the fluctuating part of the 
wind stress forcing on the top fluid layer, and the Ekman constant r 
measures the rate for vorticity decay due to the friction in the bottom 
Ekman layer. 
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1. Introduction 

The continuously stratified, three dimensional (3D) baroclinic quasigeostrophic 
flow model describes large scale geophysical fluid motions in the atmosphere 
and oceans. This model is much simpler than the primitive flow model or 
rotating Navier-Stokes flow model. When the fluid density is approximately 
constant, this model reduces to the barotropic, single-layer, two dimensional 
(2D) quasigeostrophic model. The two-layer quasigeostrophic flow model, in 
which the fluid consists of two homogeneous fluid layers of uniform but dis- 
tinct densities p\ and p2, is an intermidiate system between the single-layer 
2D barotropic flow model and the continuously stratified, 3D baroclinic flow 
model. 

The two-layer quasigeostrophic flow model has been used as a theoretical and 
numerical model to understand basic mechanisms, such as baroclinic effects 
[||], wind-drven circulation || |], the Gulf Stream f24|, fluid stability || 



and subtropical gyres [BQ, 28], in large scale geophysical flows. 



We consider the two-layer quasigeostrophic flow model (|29|1, p. 423): 

lit + J(ih , 9i + PV) = » A Vi +f + W, 

lit + J(ip2, 12 + Py) = vA 2 ip2 ~ rAip 2 , 

where potential vorticities q±(x,y,t), q2(x,y,t) for the top layer and the 
bottom layer are defined via stream functions ipi(x,y,t), ip2(x,y,t), respec- 
tively, 

Qi = AV>i - Ft ■ (V>i - ip 2 ), 

q 2 = Aip 2 -F 2 -{ip 2 -'<pl)- 

Here x, y are Cartesian coordinates in zonal (east), meridional (north) di- 
rections, respectively; (x,y) G O := (0, L) x (0, L), where L is a positive 



number; Fx,F 2 are positive constants defined by (see also p. 87) 

tp _ fo Po 



Fo 



ghi P2- P\ 
. fo Po 



gh 2 P2 - Pi 

with g the gravitational acceleration; h\,h 2 the depth of top and bottom 
layers, pi,p 2 the densities (p 2 > pi) of top and bottom layers, respectively; 
and L, po the characteristic scales for horizontal length and density of the 
flows, respectively; fo + j3y (with /o,/3 constants) is the Coriolis parameter 
and (5 is the meridional gradient of the Coriolis parameter, and v > is vis- 
cosity. Note that r = / 2 n^j^j is the Ekman constant ( [30], p. 29) which 
measures the intensity of friction at the bottom boundary layer (the so- 
called Ekman layer) or the rate for vorticity decay due to the friction in the 
Ekman layer. Here Se = \/ 2iy / fo is the Ekman layer thickness ( |]29|| , p. 188). 
Moreover, J(h, g) = h x g y — h y g x is the Jacobi operator and A = d xx + d yy is 
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the Laplacian operator. Finally, f(x,y,t) is the mean (deterministic) wind 
forcing with average zero: J fdO = 

An important part of the above equation is the white noise term W, de- 
scribing the fluctuating part of the external forcing in the top fluid layer, 
see Hasselmann [^] and Arnold (2| . An example of these short time scale 
influences is the weather variability or wind forcing. A white noise is given 
as the generalized time derivative of a Wiener process W(t) in a function 
space. 

We assume periodic boundary conditions for tp in x and y with period L. In 
addition, we impose that 



o 



ipdO = 0. 

We also assume an appropriate initial condition 



q(x,y,0) = q (x,y). 

Stochastically forced QGE has been used to investigate various phenomena 
in geophysical flows |27], 21, ||]. Recently Salmon |?2| introduced 



some generalized two-layer ocean flow models. 
It is our aim to study the long time dynamics of the stochastic differential 
equation with random coefficients but without white noise which we obtain 
after a random coordinate transformation. The solution of this differential 
equation generates a random dynamical system. The structure of this differ- 
ential equation enables us to prove the dissipativity of the system. In order 
to investigate the long time dynamics we will apply the method of determin- 
ing Junctionals in a version for random dynamical systems which is based on 
the convergence in probability, see Chueshov et. al It follows that the 
asymptotic long time behavior is determined by the asymptotic behavior of 
finitely many functionals. For instance, for these functionals we can choose 
the Fourier modes given by the eigenfunctions of the Laplacian. According 
to this property we can further show that the functionals defined on the top 
fluid layer alone can determine the asymptotical behavior of the complete 
two-layer system, when the fluid parameters and the Wiener process satisfy 
a certain condition, such as the Ekman constant r is sufficiently large and 
the trace of the covariance operator of the Wiener process is controled by 
some upper bound. Note that the generalized time derivative of the Wiener 
process models the fluctuating part of the wind stress forcing on the top 
fluid layer, and the Ekman constant r measures the rate for vorticity decay 
due to the friction in the bottom Ekman layer. 

We first recall some basic facts in random dynamical systems in Section 
2. In Section 3, we establish the well-posedness of the stochastic two-layer 
quasigeostrophic model by transforming it into a coupled system of random 
partial differential equations. The main results on asymptotic probabilis- 
tic determining functionals are presented in Sections 4 and 5. Finally we 
summarize our conclusions in Section 6. 



1 



IGOR CHUESHOV, JINQIAO DUAN, AND BJORN SCHMALFUSS 



2. Random dynamical system and determining functionals 

It is our goal to study the long time dynamics of ([l]) which is influenced by 
random forces. Appropriate tools to treat this equation are given by the 
theory of random dynamical systems. 

A random dynamical system consists of two components. The first compo- 
nent is a metric dynamical system (Q, T , P, 9) as a model for a noise, where 
(f2, T , P) is a probability space and 9 is a T ® B(M), T measurable flow: we 
have 

#o = id, t + T = Ot o T =: 6 t 6 T 

for t, t € M. The measure P is supposed to be ergodic with respect to 9. 
The second component of a random dynamical system is a i3(R + ) (g) T <8> 
13(H), ,8(£f)-measurable mapping ip satisfying the cocycle property 

ip(t + t, u, x) = (p(t, 9 T u>, (p(r, u, x)), 99(0, uj, x) = x, 

where the phase space H is a separable metric space and x is chosen arbi- 
trarily in H. We will denote this random dynamical system by symbol if. 
A standard model for such a noise 6 is the two-sided Brownian motion: Let 
U be a separable Hilbert space. We consider the probability space 

(C (R,U),B(C (R, U)),P) 

where Co(K, U) is the Frechet space of continuous functions on IR of uniform 
convergence on compact intervals which are zero at zero and B(Cq(W, U)) 
is the corresponding Borel <r-algebra. Suppose that we have a covariance 
operator Q on U. Then P denotes the Wiener measure with respect to Q. 
Note that P is ergodic with respect to the flow 

(3) 9 t u) = u(- + t) -uj(t), for u € C (K, U) 
which is called the Wiener shift. 

A main source of a random dynamical system is a random differential equa- 
tion. For example, let us consider the following evolution equation in some 
Hilbert space 

du 

(4) — = F{u,e t u;), u{0)=x, 

over some metric dynamical system P, 9). If (0) is well-posed for 

every uj G and solutions u(t,uj;x) depends measurably on (t,tu,x), then 
the operator 

if : (t, uj, x) — > u(t, u; x) 

defines a random dynamical system (cocycle) <p. For detailed presentation 
of random dynamical systems we refer to the monograph by L. Arnold 

Motivated by deterministic dynamical systems we introduce several useful 
notions from the theory of random dynamical systems. 
An uj depending closed set B contained in the separable Hilbert space H is 
called random if the mapping uj — > sup^g^^) — y\\ji is a random variable 
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for any y G H. 

A random dynamical system is called a dissipative if there exists a compact 
random set B which is forward invariant: 

ip(t,u,B(u})) C B(9 t uj), t > 0, 

and which is absorbing: for any random variable x{oj) G H there exists a 
t x (u>) > such that if t > t x (u>) 

ip(t,u>,x(uj)) G B(0tu)). 

For the following we need tempered random variable: A random variable x 
is called tempered if 

t -► \x(9 t L0)\ 
is only subexponentially growing: 

log + 

limsup [-j = (J a.s. 

This condition is not a very strong restriction because the only alternative 
is that the above lim sup is oo which describes the degenerated case of sta- 
tionarity, see Arnold [jj], page 164 f. 

Our main purpose will be to estimate the degree of freedom of the long time 
dynamics of a random dynamical systems which stem from the two-layer 
flow problem of the ocean introduced above. We will apply the theory of 
determining functionals to estimate the degree of freedom. Since dynamics 
of this model is influenced by random forces we will apply the theory of 
determining functionals related to the convergence in probability. 
We now give our basic definition: 

Definition 2.1. Let V be a Banach space which is continuously embedded 
into H . Assume that there exists r > such that ip(t, ui, x) G Lf oc (r, +oo; V) 
for almost alluj G £1 and x G H . A set C = {lj, j = 1, • ■ ■ , iV} of linear con- 
tinuous and linearly independent functionals on V is called asymptotically 
determining in probability if 
rt+l 

(P) lim / max \L((p(r, uj, xi(uj)) — (p(r, u, X2(ui)))\ 2 dr = 

t->coJ t j 

for two initial conditions x\(lo), X2(u)) G H implies 

(P) lim \\ip(t,uj,xi(uj)) - tp(t,uj,x 2 (uj))\\H = 0. 

t— >oo 



The theory of determining functionals was started with the papers [1£] 
and [^] devoted to 2D Navier-Stokes equations. Now this theory is well- 
developed for the deterministic systems (see, e.g., |lO|, 11, [l8|, 2C, 25] 



and the references therein). Some results are also available for stochastic 
systems (see, e.g., §, 0, H, U). One of the main advantages for this 
theory is the possibility to localize spatial domains or parameters which are 



6 



IGOR CHUESHOV, JINQIAO DUAN, AND BJORN SCHMALFUSS 



responsible for long-time dynamics. The existence of a finite number of de- 
termining functionals means that the long-time behaviour of the system is 
finite-dimensional. Moreover the values of these functionals on solutions can 
be intertreted as a result of some measurament of the system. Prom applied 
point of view the finiteness of number of determining functionals means that 
we need only a finite number of devices to observe completely all dynamics 
of the system. 

In this paper we rely on the following result. 



Theorem 2.2. Assume that a random dynamical system ip has an absorbing 
forward invariant random set B inV such that sup^g^^ \\x\\ v is bounded by 
a tempered random variable and t — ► sup xeB ^ tUJ ^ \\x\\y is locally integrable. 
Let C = {lj : j = 1,...,N} be a set of linear continuous and linearly inde- 
pendent functionals on V . Suppose there exist constants cc > and Cc > 
and a measurable function Ic{x\,X2,lo) which maps V x V x £1 in R such 
that Ic(xi,X2,uj) > — cc uniformly with respect to (xx,X2,lj) £ V x V x 17 
and for xi(lu), XvipS) G B(oj) we have 

V(t,u)-V(s,uj) <C C - I N c {r,uj)dT 
(5) t ^ 

+ J l C ((f(T,UJ,Xi),ip(T,UJ,X2),9 T Uj) ■ V(T,Oj)dT 

for all t > s > 0, where 

V(t,u) = \\<p(t,u),xi) - <p(t,u,x 2 )\\ 2 H , 
Mc{t,uj) = max \L((f(t, u, X\) — <p(t, lj, X2))\ 2 ■ 

j=l,...,N 



sup / Ic(<p(t,lj,xi), lp{t,lj,X2), 9 T Lj)dr > < 



Assume that 

1 I 

-E I 

for some t > 0. Then C is a set of asymptotically determining functionals 
in probability for random dynamical system ip. 

Proof. It is easy to find from relation (|5|) that 

V{t,Lj) < V(Q,Lj)eti l ^ s '^ ds + C c f M c {T,Lo)e^ lci - s ^ ds dT, 



where l(t,u) = lc{ l p{t,LO,xi),p{t,LO,x 2 ),9tLo). Therefore we can apply the 
argument given in the proof of Theorem 2.2 fl3||. □ 



In the next section, we return to the two-layer quasigeostrophic flow model. 

3. Mathematical setup and well-posedness of the random 

two-layer flows 

In the following, Lp er , Hp er for s G 1 are the standard Sobolev spaces of 
L-periodic functions with the zero mean value, i.e. J ipdO = 0. Let (-,-)o 
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and || • Ho denote the standard scalar product and norm in L 2 er , respectively. 
Every element u(x, y) G Lp er can be represented in the form 



u(x,y) 



where the Fourier coefficients Uj possess the property Uj = u-j (bar denotes 
the complex conjugation) and 

IMIo := / M 2 ^0 = \uj\ 2 < oo. 

Note that L~ 1 exp {i-£-{jxx + j 2 y)} is the eigenelement of the —A with the 

eigenvalue Ai • (jf + j 2 ), where Ai = ( 2 f-) 2 i s the smallest eigenvalue. The 
norm in FL per is defined by the formula 

\\nf s := f |(-A) f u| 2 dO = A? £ (i? + iDl^f- 
ieJPjyo 

It is clear that 

(6) l|Vit||s := H^ttll^ + ||3 y n||g = ||u||f+ g , s£R. 

We also denote Lp er = L 2 peT x I? per and Hp er = i/p er x #p er . 
We work on the phase space H~ er with the scalar product 

(9,9)* = h 1 (Vip 1 ,V^ 1 ) + h 2 (X7ip 2 , V^ 2 )o + p(^l - V , 2,^i - ^2)0, 

where 9 = (91,92), 9 = (91,92) and ^ = (ipufa) G H* er , ^ = (^1,^2)- The 
relation between (7 (resp. 9) and ^ (resp. V) is defined by (||). Here we also 
use the notation 



fo Po 



P 

9 92 ~ Pi 

Note that F\h\ = F 2 h 2 = p. The norm induced by this scalar product 

||9ll* = (9,9)* = /iillVVillo + ^HV^IIo+PllV'i -^Ho 
is equivalent to the usual norm on H" 1 ,.. Moreover, we have the estimate 

(7) MWillo + MW2II0 < Il9ll* < ao (MWi||§ + MV^||g) , 
where 

2p 2 

ao = 1 + t ^-77 — t-^ = 1 + — max{Fi, F 2 }. 

X\ mm{n,i, h 2 \ \\ 

To treat the nonlinearity we need the following lemma: 
Lemma 3.1. The operator Jacobian verifies, for u,v,w in Hp er 

(8) J(u,v) = -J(v,u), (J(u,v),v) = 0, 

(9) (J(u,v),w) = (J(v,w),u) . 
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Moreover the following estimates hold: 

(10) \(J(u,v), Au)q\ < co||A«||o • ||Vw|| ■ ||Au|| , u, v £ H 2 er ; 

(11) \(J(u,v),w) \ < ci||A«|| • ||Au|| • |H|o, u, v E H 2 er , w G L 2 er ; 

(12) \(J{u,v),w)o\ < cjUVttllo • ||Au|| • ||Vu;|| , u, w e H* er ,v e 

_ i 

//ere above cq = 2 + (v2 • vr) _1 and ci = coA x 2 . 
Proof. We start with (|Toj) . It is easy to see that 

J(it,w) • Au =- ■ v y {d x (ul - u 2 ) + 2d y (u x u y )} 

+\ • v x {d y (ul - u 2 y ) - 2d x (u x u y )} . 
Consequently relying on property @ we have 

\(J(u,v), Au) \ <^ • ||%||i (\\d x (u 2 x -u 2 y ) || _! + 2\\d y (u x u y ) ||_i) 
+\ ■ Kill [\\d y (ul - u 2 y ) ||_i + 2||da. (u x u y ) ||_i) 

<^ • (Kill + Kill) • (KIIl* + IKIli 4 + 2 K%||o) 
< (Kill + IKHi) • (IKIII4 + \\uy\\ 2 L 4) 
<V2 - \\Av\\ ■ (\\u x \\ 2 L 4 + KH^) • 
Using the inequality (see, e.g., Jig ]) 

\\u\\ L 4 < ai ■ \\u\\ 1/2 ■ ||Vn|| 1/2 with a t = ((2vr) _1 + \f2 S j* 2 , 
we obtain @. 

In a similar way we can establish the estimates 

UK \ \ 1/ IIV7 II 1/2 II A 11 1/2 MV7 1 1 1/2 I, A n 1/2 n n 

\{J(u,v),w)o\ < collvullo' • ||Ati|| ' • || Vv\\ ' • ||Av|| ' • \\w\\ , 
\(J(u,v),w) \ < co||u||o /2 • ||V-u|| ( 1 ) //2 • ||Vu||j!/ 2 • ||Av||o /2 • ||Vu>|| , 



which easily imply (11) and (|i2|). □ 

We now transform the two-layer model (||) containing white noise into a sys- 
tem of coupled random partial differential equations. The reason for taking 
such a transformation is that we need some particular a priori estimates for 
the trajectories of the system. Often these a priori estimates can be calcu- 
lated by the Gronwall lemma. For equations containing a white noise there 
exists no Gronwall lemma. But we can use this technique for the trans- 
formed random partial differential equations. 

For this purpose we introduce an Ornstein-Uhlenbeck process rj(x,y,t,uj) in 
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Lp er . This process is defined by the solution of the following linear stochastic 
partial differential equation 

(13) r] t = v(k + 1)At/ + W, 

with periodic boundary condition and with some initial condition x S F? er 
where k > is a free control parameter and W is a Wiener process in Lp er . 
We suppose that the covariance operator Q of this Wiener process has a 
finite trace. As mentioned in Section || such a Wiener process generates 
a metric dynamical system denoted by (f2, F, P, 9) where 6 is the Wiener 
shift introduced in (||). It is well known that this equation has a stationary 
solution which is generated by a Gaussian random variable rj in Hl er . In 
particular, the mapping 

t -> r/(9 t u) G Lf oc (-oo, oo;H^ er ) 

solves this equation. For moments of n we obtain: 

We introduce new variables 

(15) qi := qi - rj, q 2 , ^i := tpi + £i, ^2 := V>2 + 6, 

where stationary process n solves the problem (^) and £i and £2 are defined 
such that the elliptic equations (|2|) remain of the same form 

gi = - Fi • (^1 - tp 2 ), 

(16) _ _ 

q 2 = Aijj 2 ~ F 2 ■ (V>2 - ipi), 

The processes £1 and £2 are solutions of the linear elliptic equations 

A£i - F 1 • (& - £ 2 ) = - rj, 
( } A&-F 2 -(&-£i)=0 

and can be presented in the form 
1 



{^(-A)- 1 + A + Ft + F 2 )- x } n 



Fi +F 2 

{(_A)- 1 -(-A + F 1 + F 2 )- 1 } ?? . 



F 2 



* F x + F 2 

Thus the processes £1 and £ 2 are smoother in spatial variables than i]. In 
fact after simple calculations we have the estimates 

(18) Ui\\s+2 < His, * = 1,2, and - £ 2 || s+2 < |M| S s G R. 
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Now we obtain the coupled random partial differential equations for new 
potential vorticities qi, q 2 : 

qit+J^i ~ f 1, Qi + V + (3y) = uA 2 ^ + / 

- vA 2 ^ - v(k + 1)A?7 , 

Q2t+J (i>2 - 6, 92 + (3y) = vA 2 ^ 2 - rA4>2 

- vA 2 & + rA&. 

We treat i], £1,^2 as known processes. Using (^7|) we have 

-uA 2 ^ - v(k + l)Arj = - zaFi(A£i - A£ 2 ) - vkAr), 

-vA 2 i 2 = - vF 2 {Ai 2 - A£i). 

For convenience, we drop the tilde and rewrite the above system. Thus we 
finally get the coupled system of random partial differential equations 

qu + J(A -ti,qi+V + !3y) =vA 2 ^) X + / 

- iaFi(A£i - A£ 2 ) - vkArj 

V J O 

q2t + J(ip2-^2,q2 + f3y) =vA ip 2 - rAip 2 

- vF 2 {A& - A£i) + rA£ 2 

with 

qi =Aip 1 - Fx ■ (ip! - ip 2 ), 

q 2 =Aip 2 - F 2 ■ (1P2 - i>i), 

in the class of L-periodic functions with initial data q(x,y,0) = qo(x,y) = 
(qoi(x,y),qo2{x,y)) G H"^, where rj is the stationary solution to and 
£1 and £2 are solves (|17|) in H 2 er . 

For the rest of the paper, we work on this coupled system of random par- 
tial differential equations for stochastically forced two-layer quasigeostrophic 
fluid system. 

The coefficients of coupled system (|l9|), (|20| ) have similar properties as the 
coefficients of the corresponding deterministic twodayer quasigeostrophic 
system (cf . || , where iV-layer model with other boundary conditions is con- 
sidered). Therefore, similar to Q, using the Galerkin method and the com- 
pactness argument we can prove the following assertion on well-posedness 
of problem (|i"9l), (p0|) in the class of L-periodic functions. 

Theorem 3.2. Let go G H"^ and f G L 2 . Then for all u £ Q and for all 

T > 0, the system ( pq j has a unique solution {q(t),ip(t)} such that 

q G C([0, T\; H-i) n L 2 (0, T; L 2 er ) n L 2 OC (0, T; U l per ). 
TTie function ip associated to q by §Z(\ ) satisfies 

V G C([0,T];Hi er ) n L 2 (0,T;H 2 er ) n L 2 OC (0, T; H^). 
The solution depends continuously on the initial condition qo G H"^. 
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By the uniqueness assertion of the last Theorem the solution t — » q(t) gen- 
erates a random dynamical system 92 with continuous H~ e ^, 9 q — > w, (?) 
on the phase space H~* . 



4. DlSSIPATIVITY OF THE RANDOM DYNAMICAL SYSTEM 

Dynamical systems generated by many nonlinear parabolic differential equa- 
tions have the dissipative property which means that there exists a compact 
set absorbing the states of the system in finite time. Usually one can choose 
such a set which is also forward invariant. 

We now construct an absorbing forward invariant set for the random dy- 
namical system generated by (|l9|). This set will be a random set. 

Theorem 4.1. There exists a compact random set B(uj) C H~ e ^ such that 

<p(t, w, B{uj)) C B{O t u) for t > 0, 
(p(t,u,q(u)) C B(9 t u) for t >t (uj,q) 

where q is a random variable with values in H~ e ^.. 

We now divide the proof of this theorem in some lemmata. We start with 
the following: 



Lemma 4.2. Let q(t) be the solution of (\19l)- Then q(t) satisfies the fol- 
lowing inequality 



d 
dt 



<d ■ IMMIIo • (MVV1II0 + MW2II8) + m(0tv), 



where 



and 



m(w) = 411^)115 + 411^)113+* 



v \ Xf mmjfti, h 2 \ J 
«i = — 7 — > 

, PHhi + h 2 ) up 2 fl 1\ rh 2 , , 

d 2 =9 -3 1 I - — V 77— I + tttt- + k vhi 

uXf Xf \hi 5h 2 J I8A1 



efc=^ll/H 2 -i. 
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Proof. Multiplying the equations fll9| ) by —2hitpi, — 2/i2i/'2> respectively, and 
then adding together, we obtain 

l\\q(t)f*+Mhi\\*Ml + MAVall!)) + 2rh 2 \\VH\l 
=2^i(J(^i-^i,gi + T/ + /3y),^i)o 

+ 2h 2 {J{ip 2 - 6, 92 + A/), ^2)0 - 2/ii (/, ^1)0 + 2hivk(Arj, ^1)0 
+ 2z/p(A£i - A£ 2 , Vi)o + 2*/p(A£ 2 - Aft, ^2)0 
- 2rh 2 (A£ 2 , ^2)0- 
We have by Lemma |3.l| 

(J(ipx,q 1 + r)),ip 1 ) = 0, (J(^>2, 92), ^2)0 = 0. 

A simple calculation shows that (J(ipi, (3y), tpi)o = 0. We now estimate the 
right hand side terms. On account of (|20l) , (|Io|) , ( ^2j) , the bilinearity of J 
and (18) we can estimate 



-2 5>i(J(&, ft), V>i)o = -2 ^ J(&, Aipi),-iJji)o + 2p( J(6 - 6, ^2), V>i)o 

i=l i=l 
2 

<2c ]T ^1^^1101^^1101^^110 + 2cip||A(^i-6)||o||VVi||o||Vi/; 2 ||o 

i=l 

<~(/n||A^i||g + MA^Ug) + 6^||»;||g(/n||V^i||g + MV^||g) 

+^b^ iiv ^ ii8+ftaiiv ^ ii8) 

^l^illA^illg + ZiallA^llgJ + dollT/llg^illV^llg+^IIV^IIg). 
Similarly, due to (0) and @ 

-2^ 1 (J(C 1 ,r ? )^i)o < 2d/n||Ae 1 ||o||r 7 [|o[|A^[|o < M^M + ^I||A^||g. 

By the Cauchy-Schwarz inequality 

2 2 

-2^/ Ji (j(e i ,^) ! ^)o<2^^/ ii ||ve i |i-iiiv^iio 

i=l i=l 

< m l\t h2) H\l + £(MAfc[|g + MAiMg). 

- 2/ ll (/^l)o < 2/ tl ||/||_ 1 ||V^i|| < ^ll/ll-i + ^IIAVillg 

+ 2/1^(^,^1)0 < g^iz/iii^iig + ^iiA^ng. 
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Using (|^) once more 

2i/p(A(& - Vi)o < 2up\\A^ - A6H-2IIAV1H0 

,,9 vh] „ . „o 

H|§ + — -|AVi|g, 



^ 9zvp 2 n n2 , vh\ 



■2^(Ai;. - - i2 ). r,)„ < ^-Nl§ + ^||Ar 2 ||ii. 

5Af/i 2 



-2r/i 2 (A6,^ 2 )o < ^INo + 2r/i2||V^2||g. 

Adding all these inequalities together we obtain the conclusion. □ 

We now consider the random variable 77 defined in Section g. Recall that 
77 depends on the control parameter k. If A; is chosen large enough then 
particular moments of 77 are small. Especially we can formulate: 

Lemma 4.3. Let W be a Wiener process in L? per with finite trace of the 
covariance. Then under assumptions 

2doa troQ WdptrpQ 

i/ie random variable 



R q (oj) — I e a o Jt " a t ;no m{Q T uj)dr 
J —00 

is finite and tempered. Moreover 

(ERlf 2 < d\ ( 3a ° ^ ( 2XlV 4d ° tr ° Q ' 



2u\\J \ ao Aiz/(/c + l) 



\2 



, ~l/4 di( tr 0<2) 1/4 d 2 tr Q 
d4 ~~ G 8 x2 2(1 1 i\2 + U 4 ~\ 71 1 T\ + "3 

is an estimate for (Em 4 ) 1 / 4 fi/ie constants C§, C4 are defined in jjL&)). 



The proof of this lemma can be found in Chueshov et. al. for an 

Ornstein-Uhlenbeck process in another Hilbert space. However the argu- 
ment given there is of a general nature. 



We now construct a set satisfying (|2l|). 



Lemma 4.4. Let R(u) := aRo(tj) for some a > 1 and Rq as in Lemma 
Ejj. T/ien i/ie closed Hp^-ball -B(0, R{uj)^ ) fullfills provided conditions 
l2i) hold. 



Proof. Using Lemma 4.2 and relation (|?]) we have 
d ,, /.m, 2 , / z^Ai 



dt 1 



< (~ + do • iMMIIo) • + 
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Let qo = q(0) and p(t,u>, ||<7o||*) be the solution of one dimensional random 
affine equation 

(23) ^ + ^p = d \\r,{e t u)f p + m(e t Lo), p(0,u, \\q \\l) = \\q \\l 
at ao 

A comparison argument gives that 

y(t,u;,qo)\\l = \\q(t)\\l < P {t,u,\\q \\l). 

Here ip is the dynamical system introduced in Section |||: f(t,oj,qo) = q(t), 
where q(t) is the solution to (|l9|) with the initial data qo. Equation (|23|) has 
the stationary solution given by t — > Ro(6tUj): 

p(t,u,R {u;)) = R {9 t oj). 

This can be checked by the variation of constants formula. This solution 
is exponentially attracting which follows from the variation of constants 
formula again: 

\R (O t uj) - p(t,u, |MI*)| = \p(t,u),Ro{v)) - p(t,uj, \\q \\l)\ 

which tends to zero exponentially fast. Indeed, it follows from ( |i~4| ) and ( |22] ) 
that for a sufficient small e > 0, 



I 



do\\r](e r u)\\ 2 dT < — 1, — >0 



/o a a 

for large t > under conditions (p2|). □ 

It remains to prove the existence of a compact set B satisfying (p3|). 

Lemma 4.5. Suppose that the random variable R(lo) is defined in Lemma 
The set 



B{uj) := </?(1,6Licj,£(O,E(0_io;)2)) 
is a compact absorbing forward invariant random set. Moreover 
(24) o;^sup{||A^i||g + ||A^i||g : (gi,©)Gfl(w)} 

is a tempered random variable (ipi and ip2 o^e defined by p?Z|j ). 



Proof. Since ip(t, u, •) is completely continuous for t > (see the regularity 
assertion of Theorem |3.2| ) the sets B(u>) are compact. Since R is a random 
variable the ball -6(0, i?a ) is a random set. The continuity of (p(t, u>, •) allows 
us to conclude that B is a random set. The construction of B ensures that 
that set in absorbing and forward invariant. The temperedness of ( p4| ) can 
be proved in the same way as in [13] for the 2D Navier-Stokes equations. □ 



Hence we have shown that the two-layer quasigeostrophic flow system is 
dissipative in the probabilistic sense. In the next section we consider as- 
ymptotic probabilistic determining functionals. 
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5. DETERMINGING FUNCTIONALS 

In this section, we compare two solutions as a way to understand asymptotic 
dynamics in the probabilistic sense. Consider the difference of two solutions 
q and q, corresponding to the stream functions ip, t/j . We set 

q = (91,92) = q-q, 

corresponding to the stream function 

ip = ip — ip. 

We get random partial differential equations for the solution difference q = 
q — q from (|l9|) : 

q u =v/\ 2 ipi 

(25) - J(ih,qi+0y)- J$i,qi)- Jtyuv) + J(Sl,Ql) , 
q 2t =vA 2 ip 2 ~ rAip 2 

(26) - J(V> 2 , 92 + 0y) ~ J$2, 92) + J(&, 92) , 

Multiplying ( p5| ) by —hitpi, and (|^) by —hi^i and adding together, we get 

(27) ljt M * = ~ ^[^H^llo + M A ^Hol " rh 2 \\VMl 

+h(J(ijji - 6,91)^1)0 + h 2 {J{4>2 - 6,92)^2)0- 

We start with 

hi{J{$\ - 6,gi),'0i)o = h +h, 

where 

h = /ii(J(^i -fi,AV>i),i/>i)o, h =p{J($i - fi, ^2)^1)0- 
By © we have 

|Ii| < oofcill^i - 6ll2||A^i||o||VV'i||o- 
Using (^) we obtain 

\h\ < cipll^i-aibllV^llollV^Ho. 

Thus 

\h\ <cipA^||^i-Ci||2||A^ 2 || ||V^i||o. 



Using the inequality ab < ea 2 + -^b 2 , we get 

h(J(A - 6,91)^1)0 < - (/iillAVillo + ^IIA^IIo) 
(28) 2 

+ ^0 ( l +jFlF2 A r 2)||^_ 6 ||2||V^||2. 



16 IGOR CHUESHOV, JINQIAO DUAN, AND BJORN SCHMALFUSS 

In a similar way we have 

h 2 (J(i>2 - 6,92)^2)0 < - (/n||A^||g + MA^Ho) 

(29) /^C 2 

+ -^(1 + FiF 2 Ar 2 )ii^2 - 6iiiiiv^ 2 ||g. 



Using (||), (||) and @ we obtain 
d 



(30) (it 



I < - uihxWA^Wl + h 2 \\Aip 2 \\ 2 ] - 2rh 2 \\Vi{j 2 \\l 

+ b (M^i " eilllllVVillo + M^2 - fclllll Wallg) 



where 60 = "^(l + F1F2X1 )• This estimate is the main point in the con- 
struction of determining functionals. It is our aim to obtain an integral 
equation like (||). To this end we first consider the case when parameter 
r > is arbitrary. Let C = {lj} 1 be a set of linearly independent bounded 
linear functionals on the space H 2 er . Assume that the set C possesses the 
property 

(31) l^llx < C c max \lj&)\ +s c U\\ 2 , 

3 

for any ifi £ Hp er with some positive constant Cc and ££. We note that 
the best possible value of the parameter ec is called the completeness defect 
of the family C with respect of the pair of the spaces H per and H 2 er (see 
[10, 111]) and the smallness of the parameter ec is crucial in the subsequent 
considerations. We refer to [1C, |ll]] for general properties of completeness 
defect and for estimates of ec for several families of functionals on Sobolev 
spaces. 

From ( |3l| ) for 1/) = (y/h~iipi, \f~h 2 i\) 2 ) we have 

(MAVillo + MAV2II0) >^ (M Wi||g + \\h 2 v,p 2 \\ 2 ) 

£ c 

-Cgr max \lj(y/h^ipi, \fh 2 t\) 2 )f . 
3=1,- ,N 

for any < S < 1 with appropriate positive constant Cgc- Therefore from 
(Q and (@) we obtain 

(32) £\\q\\l<l%\j>,u)\\q\\l + C s ,c max MVhfa, Vhfa)\ 2 , 
at 3=1,— ,Jy 

where 

l c l \i),uj) = -- — y + 6 max {H^i - , ||^ 2 - 6II2} • 
aoe c 

The application of Theorem |2.2| gives the following assertion. 

Theorem 5.1. Let £ = {lj : j = 1, • • • , N} be a finite set of linearly inde- 
pendent continuous functionals on the space H 2 er . We assume that this set 
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satisfies @). Let $2%) hold. If 



^ vy/mm{hi, h 2 } 
\/2aoboY, 

where 

(33) £ = d (E\\ri\\fy (ER 2 )^ + Em + v min{/ii, /i 2 }E||V |2 



o 



anc? m ; involving the Wiener process through the Ornstein-Uhlenbeck pro- 
cess rj, is defined in Lemma ^.1 , f/ien C hl,h2 = {l^ 1,fl2 ; j = 1, • • • , AT} 

where lj 1,fl2 (i/Ji, 1P2) '■= Ijiy/hl^l, y/h2^2)i lj & C is a set of asymptotically 
determining functionals in probability for the stochastically forced two-layer 
quasigeostrophic fluid system JT^j, 



Proof. We integrate (|32|), The conditions of Theorem (|2.2| ) are satisfied if 
EZ£(i) < 0. can be estimated by 



-v- 



1 - 5 2vb 



a§e\ i/min{/ii, /12} 
On account of Lemma IO 



(7» 1 ||A^ 1 ||g + 7» 2 ||A^ 2 ||g)+26o||T;||g. 



i£B(w) * 

has the bound 



sup -v / (^i||A-0i(T,u;,a;)||g + /ti||A-0 2 (T,a;,a;)||g)dT 



\\rj(e T u)\\lR(9 T u) + m(^ T w))dr 



since 5 is forward invariant. Note that the expectation of R/t can be made 
arbitrarily small if t is large. Calculation the expectation of the last expres- 
sion and choosing t sufficiently large yields the conclusion. □ 



Remark 5.2. i) It is easily seen that there exists a set £ hl,h2 of linearly 
independent linear bounded functionals which are determining with 

respect to q. Indeed, q and ip are connected by a linear homeomorphism A 
from Lp er to ~R 2 eT defined by ( |20| ) such that we can set f 1 . 1 ' 112 = l^ 1,fl2 o A 

where C = {lj}i defines the set of determining functionals introduced in 
Theorem |5.1| . 

ii) If h\ = hi then the set C itself is determining in the sense of Definition 
P 

iii) Assumption (|22] ) holds, if 



(34) ? — < 1, k + l> — , 
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for example. In this case it is easy to see that (Ei?^) 1 ^ 2 < Ad\ (^Sl) (we 
choose a = 4/3 in Lemma 4.4). Therefore using (|l4|) we obtain 



,1/2 troQ / 2 / 3a x 2 



' \{k + l) ' [ 4d ° 4 \2^l + " min i h ^ h ^ ) +^4. 



Using ([34|) again we have the estimate 

/c + 1 V 4aoao Aiz/ 



Therefore under conditions (34) in the deterministic limit troQ — > we 
obtain estimate 



(35) e £ < — i^- 1 lL iJ • <^ (Ai + 2max{F 1 ,F 2 }) 1 + 1 " * 



6 



iCqV^i I V A i / 



-1/2 



If this estimate holds, then functionals {Ij 1 ' 2 } are determining for the deter- 
ministic two-layer model. Moreover under condition (|35j) these functionals 
are also determining for our stochastic two-layer model(|l|) and @ provided 
the noise parameter troQ is small enough. 

At the end of this section, we will discuss appropriate values for fluid pa- 



rameters. We will estimate the sufficient condition in Theorem 5.1, namely 
the inequality (|35|), in terms of fluid parameters. 

We further show that the long-time dynamics of the two-layer geophsysical 
fluid system is determined by the long time dynamics of the top layer alone , 
when the fluid parameters and the Wiener process satisfy certain conditons. 
In this case the long time dynamics will also be determined by only finitely 
many functionals. 

Let C = {lj : j = 1, N} be a finite set of linearly independent continuous 
functionals on the space H 2 er of the top stream functions. We assume that 



(36) ||VVi||o <e c - ||A^i||o + C £ - max l^i)|, Vi € H 2 

3=1, ...,N 

where Cc > is a constant depending on C As above we have 

fci||AV>i||o ^ ^-MW>lHo - Cs,c h i max |^(V>i)| 2 - 
£ c j=l,-,N 

for any < 5 < 1 with appropriate positive constant Csc- 

Thus under the condition ( |3"6| ) and using the inequality HV^Ho ^ Ai|| Allien 

the estimate (pOj) for the solution difference q — q becomes 



^IMI* + ^1^2— llV^illi + ("M + 2r)/i 2 ||VV>2||o 

< b max { H^i - , ||-i/)2 - 6111} Ikll* + C 8,C max ^(/iiV-i)) 2 . 

3=1,— ,N 
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Thus we obtain 
d 



^\\q\\l<lW$,u)\\q\\l+C s ,c max Mh^f 
ax ~ 3=1,— ,N 



where 



If (w, uj) = — mm < v T , > 

[ a e c a J 

+ feo max {ll^i - , \\tj} 2 ~ 6 111} 

Again applying Theorem |2,2| we obtain the main result. 



Theorem 5.3. Let C = {lj : j = 1,...,N} be a finite set of linearly in- 
dependent continuous junctionals on the space Hp er . We assume that this 
set satisfies the following condition involving only with the top fluid layer 
dynamical variable, i.e., the stream function ipi{x,y,t) 

(37) HWiHo <e c - ||AV>i||o + C£- max M^i)\, ^i^B 2 

3=1,... ,N 

If 

vmm{hi, /12} 

, yvAi t AT) { 

-c 

where £ is given by fiS^j), then £ is a set of asymptotically determining func- 
tional in probability for the stochastically forced two-layer quasigeostrophic 
fluid system (|T^j, Ml 



(38) E<mm^,M 1 + 2r) 



The proof of this theorem is the same as the proof of Theorem 5.1 



Relation 
(39) 

and 
(40) 



m holds if 



£ < (1/A1 + 2r) 



^min{/ii, /12} 



26 ao 

The parameter £ depends on r via c?2- Therefore it is not clear whether 



( flOP holds for some r. However as in Remark 5.2 in the deterministic limit 
ItqQ — > the estimate (4C) turns into the relation 



(41) v\ 1 + 2r > 



36cg 



1 + — max{Fi,F 2 } 
Ai 



i/ 3 Ai min{/ii, /12} 
This observation leads to the following assertion 



1 + 



F1F2 



Corollary 5.4. Assume that (41) holds. Let C = {lj : j = 1,...,N} be a 



finite set of linearly independent continuous functionals on the space H? 



per • 



We assume that this set satisfies (31_) with the parameter ec satisfying (fff^). 
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Then there exists 7 > such that C is a set of asymptotically determin- 
ing functionals in probability for the stochastically forced two-layer quasi- 
geostrophic fluid system ftll\), ( ffZ| j provided tioQ < 7. 

Consequencely, the asymptotic probabilistic dynamics of the stochastically 
forced two-layer quasigeostrophic fluid system is determined only by the top 
fluid layer. 

Remark 5.5. Note that the main task is to prove the existence of deter- 
mining functionals of the main equation (|l|) . It is easily seen by the structure 
the transformations (|l~5| ) a set C is determining in probability for (|i~9|) if and 
only if the same set C is determining in probability for ([!]) ; see |13| . 



As Theorems |5.1| and |5.3| show the problem of describing of finite families 
of determining functionals is reduced to the study of sets of functionals for 
which the estimate (|3ll) (or (|36|) ) holds with ec small enough. It is also 
important to calculate the best possible value for ec for the given family of 
functional. For this parameter there is the estimate from below depending 



only on a number of functionals [ 10 1 and this estimate coinside with ec when 
functionals are modes (see Example 1 below). We refer to [10, [ll|] for further 
discussion concerning an optimal choice of families of functionals with prop- 
erties like (EH) or (pq) . We also note that the deterministic counterparts of 
Theorems jfTI and |5.3| were proved in for other boundary conditions. 

Example 1 (Determining modes) Let {ei(x,y)}^ =1 be the basis of eigen- 
functions of —A with the periodic boundary conditions in O such that the 
corresponding eigenvalues possesses the property 

(42) < Ai < A 2 < . . . . 

We note that every eigenfunction has the form L _1 exp{^(j'ix + 32V)}- 



However we numerate them such that (f42|) holds. Suppose Cn = {lj 
j = 1, N} is the set of the functionals on Hp er of the form 

lj(u)= / u(x,y)ej(x,y)dO, j = l,...,N. 
Jo 



Then one can prove (see, e.g. pi) that estimate (36) holds for Cn with 



— 1/2 

££jv = ^n+i ana - * ms value is the best possible among all families of func- 
tionals consisting of elements. Since A;v ~ cqNL~ 2 for A^ large enough, 
we have that ec N ~ cqL/^/N with some absolute constant cq. Therefore 
we can use Theorem [5.3| to estimate number of determining modes. We 
also note that using the family Cn of the functionals on H^ er we can easily 
construct a family Cjy of the functionals on H5; er = Hp er x Hp er such that 
(|3l]) holds with ec N = ^{i- 

Example 2 (Determining nodes) Let us consider the nodes 
(xi, yj ) = -^=-(i,j)eO, l<i,j<VN, 
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and define functional kj on H^ er as 5-functions at (xi,yj), i.e by the for- 



mulas kj(u) = u(xi,t/j). Let C = {kj}- One can prove (see, e.g. [|To|] or 
p5|) that estimate (B^) holds for C with sc = cqL/^/N where cq = \ is an 
absolute constant |25| j. Thus S£ has the same order for large N as is the 
best possible value £c N for families of functionals consisting of N elements. 

We also note that the estimates for the completeness defect ec and the 



Ekman constant r given in Theorems 5.1 and 5.3 are rather crude. In 
the two-layer quasigeostrophic model considered above, the parameter v 
is the molecular viscosity. For ocean water, v = 10 _6 m 2 s _1 . However, 
fluid turbulence at small scales can act as an extra dissipative mechanism, 
thus calling for the substitution of the molecular viscosity by a much larger 
eddy viscosity. For example, in the two-layer quasigeostrophic model flow 
simulation in [p8| , the eddy viscosity (we still use the same notation as the 
molecular viscosity) is taken as v = 50m 2 s . The gravitational acceleration 
g = 9.81ms -2 . 

At mid-latitude (45 degreeiV), f = 8x lCT^ 1 and (5 = 2.3xl0 _u m- 1 s- 1 . 
For large-scale flows at mid-latitude, such as the Gulf Stream in the Atlantic 
ocean, the horizontal spatial scale L is at the order of 1000A;m. Moreover, 
layer depth h\ and hi are at the order of 500m each for large scale flows 
such as the Gulf Stream at the mid-latitude in the Atlantic ocean. For the 
eddy viscosity 50m 2 s _1 in pH , the Ekman constant r is then at the order 
of 10 _5 s _1 . The number r is large when the eddy viscosity is taken to be 
large. 

The ocean water mean density po is about 1025/cgm~ 3 or 1.025<7cm _3 . pi 
and p2 are in the vicinity of po. The density difference p2 — pi ~ 25/cgm~ 3 . 
But it is this small density difference that in turn determines pressure dif- 



ferences and thereby drive the ocean circulation [34]. 

The mean wind forcing /, i.e., the deterministic part of the curl of the wind 
stree on the top fluid layer is usually taken as a stationary, or being time- 
averaged and even also zonally averaged, sinusoidal function. For example 



2ttt . 27ry 
PohiL L 

where the wind tension tq is of order 1 dyne/ cm? or of order 0.1 N/m 2 . 
With this mean wind forcing and physical parameters specified above, the 



inequality (35) turns into the estimate ec < 1-35 • 10 _3 m and therefore in 
the case of Examples 1 or 2 for the number N of functionals we obtain the 
estimate N ~ 10 18 . Thus Theorem 5.1 should be only considered as an 



qualitative assertion about finite-dimensionality of the long-time behaviour 
of the stochastically forced two-layer quasigeostrophic fluid system. As for 
Theorem |5.3| , a similar calculations show that the condition ( |4l| ) can be 
valid under some special choice of parameters. Therefore Theorem |5.3j only 



22 IGOR CHUESHOV, JINQIAO DUAN, AND BJORN SCHMALFUSS 

predicts the possibility of a situation when the bottom fluid layer is slaved 
by the top layer. 



6. Summary 

We have considered asymptotic probabilistic dynamics of the stochastically 
forced two-layer quasigeostrophic fluid system. We first transformed this 
system into a coupled system of random partial differential equations and 
then show that the asymptotic probabilistic dynamics of this system depends 
only on the top fluid layer, provided that the Wiener process and the fluid 
parameters satisfy a certain condition, i.e., the inequality (|38|); see Theorem 



5.3, In particular, this latter condition is satisfied when the trace of the 
covariance operator of the Wiener process is controled by a certain upper 
bound (see Corollar y |5.4| ) and the Ekman constant r is sufficiently large 
(see the inequality (|4lD). Note that the generalized time derivative of the 
Wiener process models the fluctuating part of the wind stress forcing on the 
top fluid layer, and the Ekman constant r measures the rate for vorticity 
decay due to the friction in the bottom Ekman layer. 
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